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LTX:  531.55:531.2 

MINIMUM-FL5EL  OPTIMAL  MLXjTIPULSE  .APSIDAL  INTERORBITAL  FLIGHTS 

907F0286A  Leningrad  ATISTNIK  LENINGRADSKOGO  L'NIV’ERSITETA:  MATEMATIKA, 
MEKHANIKA,  ^ASTRONOMIYA  in  Russian  No  1,  Jan  90  {manuscript  received 
6  Apr  89)  pp  63-67 

[Article  by  S.  N.  Kirpichnikovl 

[Textl  The  problem  of  finding  minimum-fuel  optimal  pulsed  flight 
trajectories  between  Keplerian  orbits  in  a  central  Newtonian 
gravitational  field  is  studied.  The  time  of  motion  along  the  orbits  is 
free.  Unlike  the  similar  problem  of  interorbital  flight,  in  which  the 
orbit  of  the  flight  from  the  initial  to  the  final  orbit  must  be 
constructed,  one  must  construct,  a  flight  trajectory,  tegun  on  the 
initial  orbit  and  which  only  intersects  the  final  orbit,  in  the 
consider€?d  problem  of  interorbital  flight,  so  that  the  equalizing  pulse 
rate  at  the  finish  point  is  not  achieved.  The  characteristic 
maneuvering  speed  is  minimized  with  the  following  assumptions.  1) 
"fragmentation"  of  the  pulses  is  not  permitted.  2)  the  total  number  of 
pulses  is  n  'n  N,  where  N  is  given  (N=  1,  2,  ...),  and  3)  the  flow  lates 
of  gases  during  all  pulses  are  identical. 

Let  us  first  consider  the  case  of  coplanar  boundary  orbits.  Using 
argument  [11,  it  is  easy  to  prove  the  following  statement,  which  permits 
one  to  limit  oneself  to  coplanar  postulation  of  the  problem. 

Theorem  1 .  The  trajectory  of  an  minimum-fuel  pulsed  flight  between 
coplanar  boundary  orbits  lies  in  the  plane  of  these  orbits. 

The  author  [2-4]  studied  the  coplanar  problem  of  finding  minimum-fuel 
optimal  single-pulse  (N  =  1)  flight.  Using  the  methods  of  these 
articles,  it  is  easy  to  prove  the  following  statements,  somevvbat  more 
general  compared  to  the  i,ndicated  articles ,  about  the  type  of  steady  and 
optimal  orbits  of  pulsed  flights  in  the  coplanar  problem. 

Theorem  2.  Any  fixed  orbit  in  the  problem  of  an  minimum-fuel  optimal 
single-pulse  flight  always  touches  the  final  orbit  at  the  finish  point. 
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Corollary.  A  limited  fixed  multipulse  orbit  in  the  problem  of  an 
minimum-fuel  ox>timal  interorbital  flight  always  touches  the  final  orbit 
at  the  finish  point. 

The  following  four  theorems  are  related  to  fixed  and  optimal 
single-pulse  (N  “  1)  orbits  in  the  coplanar  problem  of  finding* 
minimum-fuel  optimal  trajectories  of  interorbital  flights. 

Theorem  3.  The  radial  pulse  may  not  be  fixed  in  any  boundary  orbits  at 
all. 

Theorem  4.  The  tangential  pulse  can  be  fixed  only  if  the  lines  of  the 
apsides  of  the  boundary  orbits  coincide  (specifically,  if  at  least  one 
of  them  is  circular),  and  single-pulse  flight  is  apsidal,  i.e.,  both 
apsides  of  the  Keplerian  orbit  of  the  flight  coincide  with  the  apsides 
of  the  boundary  orbits. 

Theorem  5.  Let  the  boundary  orbits  have  a  common  line  of  apsides  or  at 
least  let  one  of  them  be  circular.  The  apsidal  maneuvers  and  only  these 
maneuvers  are  fixexi  single-pulse  maneuvers  in  the  considered  extreme 
problem. 

Theorem  6.  Let  the  difference  of  the  angular  distances  of  the 
peri centers  of  the  boundary  orbits  be  a  multiple  of  n  or  at  least  one  of 
them  is  circular.  If  the  final  orbit  is  located  entirely  outside 
(inside)  the  initial  orbit,  the  global  optimal  maneuver  of  a 
single-pulse  flight  is  an  apsidal  maneuver,  which  ends  at  the  pericenter 
(apocenter)  of  the  final  orbit. 

A  system  of  necessary  conditions,  derived  in  [2],  is  used  in  the  proof 
of  these  theorems  and  the  case  of  intersecting  boundary  orbits,  trivial 
in  the  problem  of  an  interorbital  flight  with  free  time,  is  eliminated 
from  the  consideration  here  and  below^  The  following  easily  proved 
theorem  guarantees  the  absence  of  abnormal  solutions  (gradients  of  the 
left  sides  of  the  coupling  equations  in  the  system  of  necessary 
conditions  are  Ji nearly  dependent)  in  the  coplanar  problem  of  an 
minimum-fuel  optimal  pulse  flight  between  nonintersecting  boundary 
orbits . 

Theorem  7 .  Abnormal  solutions  are  possible  in  any  extreme  problem  of 
constructing  interorbital  single-pulse  flights  or  two-piilse  flights  with 
free  time  only  upon  touching  the  toundary  orbits.  These  solutions  are 
characterized  by  touching  all  three  orbits  (the  initial,  intermediate 
and  final  at  a  common  point — the  point  of  the  start  and  finish). 

It  was  established  above  for  coplanar  boundary  orbits  with  coincident 
apsidal  lines  that  fixed  and  globally  optimal  single-pulse  flights  are 
ax^sidal  .  Let  us  study  in  this  regard  fox'  the  boundary  orbits  below  the 
pi'ob.lem  of  the  minimum- fuel  optimality  of  apsidal  multipulse 
i nterorbi tal  f 1 ights . 
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Let  us  continue  as  before  to  assume  that  the  problem  is  coplanar.  It  is 
easy  to  see  that  trajectories  u'ith  apsidal  tangential  connections  of  all 
Keplerian  branches  and  the  maximum  cases  of  these  connections — apsidal 
tangential  drifts  to  infinity  and  returns  from  infinity  along 
fjarabolas — satisfy  the  system  of  necessary  conditions  in  the  case  ol 
boundary  orbits  with  common  line  of  apsides  of  the  corresponding 
considered  extreme  problem.  Let  us  also  call  the  trajectories  apsidal 
due  to  the  brevity  of  them.  The  following  statement  determines  the 
sufficient  conditions  of  global  optimality  of  apsidal  trajectories. 

Theorem  8.  If  the  orientation  of  one  of  the  boundary  orbits  is 
arbitrary  in  the  coplanar  problem  of  interorbital  pulsed  flight  (for 
example,  if  at  least  one  of  them  is  circular’),  the  character ist 3 c 
maneuvering  speed  assumes  tlie  least  value  on  the  apsidal  trajectories 
when  the  lines  of  the  apsides  of  boundary  orbits  coincide.  If  the  final 
orbit  is  lo<rated  outside  (Inside)  the  initial  orbit,  (he  globally 
optimal  maneu^-er  with  limit.ed  flight  trajectory  terminates  at  the  ^ 
pericenter  (apocenter)  of  the  final  orbit,  wbile  the  difference  of  the 
angular  distances  of  tlie  f)ericenters  of  the  boirndary  orbits  should  be 
equal  to  (2k  +  1)7,  Vvhere  k  is  an  integer,  for  single-pulsed  flight. 

Common  geometric  proirerties  of  transversally  [51  intersecting  curves  on 
a  plane  and  also  some  results  of  [1,  6-8]  and  methods  used  there  to 
study  the  problems  of  interorbital  flights  are  usfid  in  proof  of  this  and 
all  subsequent  theorems. 

Thus,  let  the  boundai-y  orbits  now  be  coplanar  and  have  a  common  line  of 
apsides.  Let  us  study  the  apsidal  trajectories  with  N  =  2.  Tivo  new 
types  are  pxDssible  besides  the  ab<3ve  apsidal  single-pulse  trajectories. 
These  are  primarily  bielliptic  apsidal  trajectories.  The  start  is  made 
at  the  apsida  of  the  initial  orbit  (polar  radius  r;)  and  the  finish  is 
made  at  the  apsides  of  the  final  orbit  (polar  radius  ri. ) ,  and  both 
corresponding  apsidal  points  lie  on  the  same  semiaxis  of  the  common  line 
of  apsides  of  the  boundary  orbits.  The  additional  (second)  .apsidal 
tangential  pulse  is  made  on  the  opposite  semiaxis  of  this  line  at  a 
point  with  the  polar  radius  p  >  0.  Finally,  both  seraiellipses  of  the 
bielliptical  apsidal  trajectory  change  to  iiarabolas  in  the  final  case  at 
p  (f>,  and  the  additional  pulse  becomes  zero  at  infinity.  The  flight 
and  its  trajectory  in  this  limiting  case  is  called  biiiarabolic .  We  note 
that  the  trajectory  of  bi parabolic  flight  with  such  a  small  increase  of 
the  characteristic  velocity  can  be  approximated  in  practice  by  a 
bielliptical  apsidal  trajectory  Vvith  rather  large  p. 

Theorem  9.  Let  the  boundary  orbits  have  a  common  line  of  apsides.  Only 
single-pulse  and  biijarabf.-)lic  flights  can  be  optimal  in  the  class  of 
flights  with  single-  and  two-pulse  trajectories  (N  =  2).  It  is 
necessary  that  n  >  Iro  for  optimality  of  a  biparabolic  flight. 

Let  us  fonnulate  the  auxiliary  obvious  statement  in  the  form  of  a  lemma. 
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Lemma.  A  pulse  drift,  maneuver,  minimal  according  to  the  characteristic 
speed,  with  given  initial  Kepleriam  orbit  to  infinity  is  al^vays  single¬ 
pulse,  and  the  pulse  is  tangential  and  is  made  at  the  pericenter  of  the 
initial  orbit. 

Let  us  now  remove  the  constraint  N  =  2.  The  following  theorem  is  easi Im¬ 
proved  bm'  the  method  of  mathematical  induction  with  regard  to  theorem  9 
and  the  lemma. 

Theorem  10.  Let  the  boundary  orbits  have  a  common  line  of  apsides. 
Further,  let  the  apsides  of  the  initial  and  final  orbits,  l.ving  on  the 
same  semiaxis  OL  of  this  line,  have  polar  radii  rj  and  rv,  respectively. 
If  maxlri/ro}  <  4  upon  two  possible  selections  of  semiaxis  OL,  some 
single-pulse  flight  is  globally  optimal  in  the  class  of  minimiun-fuel 
optimal  flights  with  pulsed  apsidal  trajectories,  and  either  biparabolic 
flight  with  start  at  the  i)ericent.er  of  the  initial  orbit  or  a 
single-pulse  flight  is  globally  optimal  in  the  opposite  case. 

Sped fi cal Ij",  if  the  initial  orbit  is  located  inside  the  final  orbit  or 
^  distance  of  the  apocenter  of  the  initial  orbit 

is  the  distance  of  the  pericenter  of  the  final  orbit),  onl.y 
single-pulse  flights  can  be  globally  optimal. 

Let  us  now  take  the  assumption  of  the  coplanarity  of  the  problem. 

Theorem  11.  If  the  biparabolic  flight  is  globally  optimal  in  the 
coplanar  problem  of  a  minimum-fuel  optimal  pulse  interorbital  flight 
with  arbitrarm"  orientation  of  one  of  the  boundary  orbits,  biparabolic 
flight  in  the  plane  of  the  initial  orbit  will  also  be  globallm’’  optimal 
in  any  problem  of  minimum-fuel  optimal  jjulsed  flight  with  arbitrarily 
give.n  arrangement  of  boundary  orbits  in  space  of  the  same  dimensions  and 
shape . 

The  statements  of  theorems  10  and  11  permit  considerable  amplification 
for  the  case  of  circular  boundary  orbits. 

Theorem  12.  A  Hohmann  single-pulse  flight  (see  f21)  is  globally  optimal 
in  the  class  of  minimum-fuel  pulsed  optimal  flights  between  circular 
both  coplanar  and  noncoplanar  l)oundary  orbits  of  radii  n  and  r^ ,  if 
ri/r-?  <  4.828  ...,  and  biparabolic  flight  is  globally  optimal  if  ri/ri>  > 
4.828....  Any  two-pulse  apsidal  bielliptic  orbit  with  /?  >  ri  is  more 
optimal  than  a  Hohmann  single-pulse  flight  for  ri/ri  >  5.879.... 

Finally,  there  is  a  critical  value  p ^  >  ri  for  any  I’j  ,  rj ,  ri/r2  € 

(4.828,..,  5.879...),  beginning  with  which  a  bielliptical  apsidal 
two-pulse  flight  is  more  optimal  than  a  Hohmann  single-pulse  flight  for 
any  p  >  and  the  economy  of  the  characteristic  speed  is  greater,  the 

greater  p ,  p  >  .  All  the  considered  maneuvers  are  made  in  the  plane 

of  the  initial  orbit. 
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Comment.  The  precise  value  of  constant  4.828...  is  equal  to 

()/2  -  1)'-  -  1.  Constant  5.879...  is  the  streatest  real  root  of  the 
equation 


—  9x^  +  I5x^  +  \9x  +  4=0. 


We  note  in  conclusion  that  constants  4.828...  and  5.8(9...  under 
conditions  of  theorem  12  play  the  same  role  as  constants  11.938...  and 
15.581...  in  the  theory  of  minimum-fuel  optimal  pulsed  interorbital 
flights  (see,  for  example,  [9],  section  3.2),  but,  unlike  the  latter, 
the  conditions  of  theorem  12  are  considerably  nons^nmmeti  ical  in  the 
considered  theory"  with  respec't.  to  radi  i  r  i  and  rv  of  the  initial  and 
final  orbits. 
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AVIATION-  AND  SPACE  TECHNOLOGY 


LHC  531.36:521.1 


ON  ONE  FORM  OF  DIFFERENTI.AL  EQUATIONS  OF  PERTURBED  MOTION  OF  GRAVITI'- 
ORIENTED  SOLID 

907F0286B  Leningrad  \ESTNIK  LENINGRADSKOGO  L'NIVEItSITETA:  M4TEMATIK4, 
MEFvHANIK4,  ,4STRONCMIYA  in  Russian  No  1,  Jan  90  (manuscript  received 
8  Jun  89)  pp  71-75 

(Article  by  A.  A.  Tikhonov] 

[Text]  Difficulties  frequently  arise  in  construction  of  the  transformed 
system  of  differential  equations  of  motion  when  solving  nonlinear 
problems  of  mechanics  using  approximate  analytical  methods,  which 
provide  transition  to  new-  variables.  Hiese  difficulties  can  sometimes 
l)e  o^•ercome  to  a  significant  degree  by  using  the  method  of  variation  of 
canonical  random  constants. 

Tills  method  is  a  convenient  means  of  constructing  differential  equations 
of  the  perturlied  motion  of  mechanical  systems,  the  state  of  which  can  be 
described  by  canonical  equations.  Thus,  the  method  of  variation  of 
canonical  random  constants  in  combination  with  averaging  the 
differential  equations  was  used  successfully  in  (1,  2]  to  study 
quasl-linear  vibrations  of  a  solid  ex-posed  to  gravitational  forces  in 
the  neighborhood  of  the  jxjsition  of  stable  equilibrium. 

A  form  of  notation  of  differentia]  equations  of  perturbed  motion  of  a 
gr-avity-oi'iented  solid,  convenient  for  study  and  which  generalizes  the 
corresponding  equations  used  in  [1,  2]  for  the  case  of  potential 
perturbing  forces,  is  found  in  this  paper.  The  rotational  moment  of  the 
solid  with  resjiect  to  its  center  of  mass,  moving  in  a  circular  Keplerian 
orbit  in  a  central  gravitational  field,  is  studied  on  the  assumption 
that  nonjXDtential  forces  that  allow  quadratic  approximation  are  acting 
on  the  body  in  a  central  gravitational  field. 

Two  coordinate  systems  wdth  common  origin  at  the  center  of  mass  of  the 
body  are  introduced  into  consideration:  c(7}( — an  orbital  coordinate 
system  (OSK)  (axis  (  is  directed  along  the  radius  vector  of  the  center 
of  mass  of  the  Ixxly  with  respect  to  the  attracting  center,  axis  (  is 
directed  along  the  tangent  to  the  orbit  in  the  direction  of  motion,  and 
axis  Tj  is  directed  along  the  normal  to  the  orbital  plane,  cxyz  is 
rigidly  bound  to  the  body  so  that  its  axes  are  directed  along  the  main 
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j. 


central  axis  of  inertia  of  the  bcxiy,  -b  k  are  the  unit 

vectors  of  axes  (.  //,  and  (  and  x,  y,  z,  respectively,  which  satisfy  the 
equalities  Sq  =  (Xii  +  aj  “K  *3^^,  ’J*  ~  =  Ti^  "f"  T2j  T3^^>  •  where 

?i.  Ti  (i=li  3)  are  direction  cosines. 

The  follo^^dng  notations  are  adopted  in  the  paper:  u'  =  e'o  +  is  the 
absolute  angular  velocity  of  the  body,  =  a'a ?/c  is  the  angular  velocity 
of  the  OSK,  =  pi  +  qj  +  rk  is  the  angular  velocity  of  the  body  with 
respect  to  the  OSK,  and  A,  B,  and  C  are  the  principal  central  moments  of 
inertia  of  the  body. 

Differential  Euler  equations  of  rotational  motion  of  the  body  diie  to  the 
effect  of  gravitational  moment  and  of  a  general  perturbing 

nonpotential  moment  Mb  ( |  Mb  |  |  Mr  | )  have  the  form 


(Aw) -j- ft)  X  (Aw)  =  Mr  +  Mb  , 


(1) 


where  A=:diag  (A,  B,  C),  Mr  =  3o)2^o  X  (A^o),  Mb  is  the  analytical  vector 
function  of  variables  a.-,  p<,  Yw  P-  ''  (i=X~3),  generally  non-autonomous . 
Equation  (1)  together  with  kinematic  Poisson  equations 


^  =  So  X  w  —  (Uo^o.  —  ’lo  X  w. 


^  =  Co  X  w  + 


(2) 


form  a  close(i  differential  system. 

If  there  are  no  j^er turbine*  moment  and  fulfillment  of  the  conditions  of 
gravitational  orientation  B  >  A  >  C  [2],  the3  body  has  a  stable  position 
of  equilibrium  Q\  =  -  73  ::  1  in  the  OSK,  taken  as  unperturbed  motion. 

It  is  determined  by  equalities  f  0  in  ’’aircraft  angles  of 

orientation  (see  figure). 

Posing  the  problem  of  finding  a  system  of  differential  equations  of 
perturbed  motion  in  a  form  convenient  for  study  of  quasi-linear 
vibrations  of  a  gravity-oriented  solid,  let  us  consider  the  nonlinear 
component  of  the  gravitational  moment  as  jierturbation  along  with 
perturbing  moment 


T 


Following  fl,  21,  let  us  construct  a  i^erturbed  differential  system  using 
canonical  variables.  Let  us  introduce  the  following  canonical 

variables :  ?.  6,  p„  p,  p^^  where  p  i  =  ^  sin  6  +  sin  9  cos  6  + 

+  COS  ®  cos  0,  />,  =  .^  =  /?3  =  -^+50)0  =  Bu>y  cos  9  —  C(o^  gin  ^  4. 

+  5u,o;  L  =  1  {^u,2  +Ba.2  +C«>2)  _  |a,2[(A-C)TH(5-C)Tl]  ^ 

Lagrange  function,  in  which  ii'x,  i’y,  i'z*  ?'i  1  and  To  are  expressed  by 
't*.  ?,  0,  >i»,  9,  0  basis  of  kinematic  formulas:  p  _ <j)  sin  0, 

^  =  ecos  9+'|isin9  cos  0,  r=  <}»  cos  9COS  0  — 0  sin  9,  91  =  — sin  0,  92  =  sin  9  cos  0. 

For  the  projections  of  relative  angular  velocity  ,  we  will  have 


V  ""  ,  «r  ““  ^  ?  +  (Pi  +P2  sin  0)  X 

X  sin  9  cos  6]  —  wg  (sin  <b  sin  9  sin  0  +  cos  <{1  cos  9),  r  =  C"*  f—  (n  ^Rio  )  v 
X  .in  ,  +  to  +p,  „  e.,  ,  ,|  _  ^  ^  ,  I  <>.+S»,)  X 


{6) 


Turning  to  the  nev:  independent  variable  u  £i/oti  let  us  introduce  the 
dimensionless  variables 


(Cu)„)-1;7,  ’)  9,  9'  =  (AiOg)-^  po  —  6'  =  (5u,o)-i  ^3. 


(1) 


Formulas  (4),  \\'hich  coincide  with  linearized  expressions  (3),  divided  by 
ii'o,  are  taken  as  the  transition  formulas  from  generalized  pulses  pi  to 

new  variables  ij)',  q)',  9'  of  the  same  order  of  smallness.  Expressing 
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a);c,  (Oy,  az,  a/,  p,,  yi  {i=l,  3)  by  q)'.  6'  and  retaining  terms  no  higher  than 

second  order  of  smallness  with  respect,  to  the  introduced  dimensionless 
variables,  let  us  transform  system  (1)  and  (2)  to  the  form 

(?')'  =  ^1,  -J*'  =  f  +  +  (1  -  Se-')  s''?  +  -  1 )  l-e. 

£8-i<p('i7  +  (0y=F2,  <p'  =  f'  +  fe, 

(f )'  +  (?)'  -  8e-'<P  (0')'  =  ^3^  6'  =  0'  +  (so-t  -  1 )  (f  <p 

where  ^  g-i  J4  ,  _  gj  _  g  _  3  ^  +  (e  —  8)  ((8  +  e-1)  y6'-eZ'B')— 

_(1  _8+e)  0^'+(8-l) <{.0]+yW, ,  /"2=8-'  [3  (s-1 )  +  (1  _  8  +  5) <??'+ 

+  (1  - s)  (<!^?  -  -!>?)]  +  3,  3  =  e-i  [(1  -  8  +  _ (8  _  1) (.1.  +  3  96)  + 

+  (e  —  1)  6f  —  8,^6']  +  M^,  Mj,  =  My  =  (8A(u^)-‘  M„y.  M^  = 

=  (e^<uo)"‘ iWszare  expressions,  calculated  with  regard  to  quadratic  terms 

and  l  =  BA-\  z^CA'K 

System  (5)  at  =  0  permits  a  canonical  form  of  notation  with 
Hamiltonian  H  =  Ho  +  H|  ,  where  HQ  =  -^{Ay''' 

+  3(^-C)0**+(i3-A)tj."]  =  2{A  —  C)l'^B^ 

+  2>Ui'9'94-S0'(1)2]  4.  (C _ 5)i{»(p0  components  of  second  and  third  orders 

of  smallness,  respectively.  Linear  approximation  of  this  system, 
corresponding  to  Hamiltonian  Ho,  i.e.,  an  unpfjrturbed  differential 
system,  has  a  general  solution  of  type 


(5) 


2 

/=1 


•i*  =  2 V" sin  tf,  9  =  2  «2/K^cos  tj,  6  =  Y'l  kY  sin  tg. 


(6) 


=  2  ®°s  "f,  ~  ~  2  —  y  2  nns  tg,  { 7  j 

‘■“1  (-1 


where  ^/(«  +  Vj)  (/=='l,  3);  Xi,  Vj  (/=1,  3)  are  canonical  arbitrary 

constants  and  kj  (,j  =  1 ,  3 )  are  unperturbed  natural  frequencies  of  the 
body: 
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Ai,2  ■= 
b  —  ^ 

Ai  =  - 


■■  ]/" a/2  +  y"a'/4  —  b,  ^3  =  1/3 (1  —  e)o  *,  a  =  1  +  3(8  —  e)  +  &/4, 
5-1  (5  _  1)  (8  -  s);  = V2  t-ikT'  = -  y  y  ^rs 

(-l)‘Se(A2_*2^  ft(l_8  +  s) 


The  inequalities  y!i.>o  (j=l,  2),  xi>0,  »X2  <  0,  aii>0,  ai2<0,  O2i<0,  a22<0. 
are  valid  in  the  region  B  >  A  >  C. 


Solution  of  system  (5)  by  the  method  of  variation  of  arbitrary  constants 

xi,  Vi  (i  3  1,  3)  is  sought  in  the  form  of  (6),  where  xi  and  vi  are  new 
unknown  functions  subject  to  determination. 

Observation.  It  follows  from  the  kinematic  equations  of  system  (.5), 
written  in  variables  xi  and  Vj ,  that  the  values  of  x'i  and  v'l  have  the 

order  0{  y/\'^ ;  )  and  0{>/x7)i  respectively.  Taking  these  estimates  into 
account  and  also  taking  into  account  that  the  differential  equations  of 
perturl>ed  motion  with  respect  to  conditions  of  postulation  of  the 
problem  should  contain  terms  no  higher  than  second  order  of  smallness 

with  respect  to  6,  "5',  ’0'  in  derivatives  (J')' ?.  (f'V 

contained  on  the  left  sides  of  the  dynamic  equations  of  system  (5),  sind 
one  should  according.] lease  only  the  follovsdng  terms  of  second  order  of 
2  2 

smallness:  cp  X  sin  — ■  0^ cos  (pk^y"2  sin 

*  frt  i^i 

uliere  and  are  determined  by  formulas  (6). 

Turning  in  system  (5)  to  variables  xi  and  vi  according  to  transform  (6) 
and  (7)  and  solving  it  with  respect  to  the  derivatives,  we  find  the 
desired  system  of  differential  equations  of  perturbed  motion  of  the 
body,  wliich  can  be  la-itten  in  the  form 


- ^  +  2(—  ly  — 02,3^102  cos 

(*/]/•*;)■*(—  l)‘['Jl.3-/3r‘AI,COST^  — 02,  3^ia2^M^sm  tj], 


•^3  =  —  ^  +  y cos 


^3.  '£'3  = 


_  SHi 


dxo 


My  sin  ^3. 


k,V‘2x, 


Here  i — 1,  2;  — kx<3^2^2i  —  ^2'^u^22»  ^2~^t'^ii^22 — ^2‘^i2®2u  and  the  expression 

for  in  variables  xf,  vi  has  the  form 
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—  (AjJCj *1*  hiX^  Xg  cos  "tj  h^Xi  |/^ x^  cos  (2 X3)  -4-  h^Xi]^ X 
X  cos  (2  ti  —  X3)  +  AjATj  Vx'i  cos  (2  xj  4-^^  h^x.  Y X3  cos  (2 Xj  —  X3)  4- 
4“  XiXnX^  cos  (Xj  4"  ''3  “  ”^2)  “f"  ^sV X1X2X3  COS  (X2  4"  "^3  ^l)  ~t" 

4"  fioY'xiX^x^  COS  (xj  X2  “■  X3)  4"  ^10  X1X2X3  COS  (xj  4"  ^2  ^3)* 

where  =  Y  Bj2  [(e8~'  —  1 )  (AaUu  —  ®2a)  32a  4“  3i«/2]  (*  =  1 ,  2) , 

h^^Ym  [(e3-‘-l)32,(ife,3„-a2,4-(-  1)P  A3“‘3i,)-aV24-(-  1)P5-‘A^X 
X  [(1  —  s)3u  — ^i32<]]  (i=l  AJin  P  =  3,  4;  i  =  2  for  P  =  5,  6), 


.  h^=Y  BIS 


2 

(s8"^  —  1 )  33  —  (  —  1  2( ^  )‘^3- 3- /[( 1  ~®)  ®  U~f^l’^2i\  + 

/=1 


4“  ( — ®)  (^12^21  —  3jj<J22)  4"  3u3j2]  (T  “  ^)’ 

h,=y  Bjb  - 1)  ®3  +  (  — [(^1*2  +  5  —  ®)  (®11®22  +  ®l2®2l)*” 

—  (1 — s)  (^1 4"  ^2)  ^11^12!  —  ®ii®i2]  (^=9.  19)> 

33  =  A1IJ11O22  4"  ^2®I2®21  2o21322" 


The  convenience  of  using  the  derived  differential  equations  of  perturbed 
motion  in  the  form  of  (8)  for  study  of  quasi-linear  vibrations  of  a 
solid  is  determined  by  the  following  circumstances;  1)  the  structure  of 
these  equations  is  such  that  their  right  sides  are  divided  into  two 
groups  of  terms,  from  which  the  first,  written  in  canonical  form, 
corresponds  to  "canonical"  perturbations,  caused  in  the  given  case  of 
nonlinear  quadratic  terms  at  the  gravitational  moment,  wbile  the  second 
group  corresponds  to  "noncanonical"  perturbations  due  to  moment 

generally  being  of  an  arbitrary  mture;  and  2)  system  (8)  has  a  form 
that  permits  direct  use  wdthout  laborious  computations  of  asymptotic 
methods  of  nonlinear  mechanics  for  study  of  it. 

Equations  (8)  are  canonical  in  the  absence  of  a  perturbing  moment  and 
coincide  with  the  widely  used  equations  in  canonical  variations  [2]. 
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[Article  by  N.N.  Berezkin,  engineer,  Yuzhtekhenergo;  and  S.S.  Kirilyuk,  C.Chem.S.,  Lvov  State  Medical 
Institute]  ^  . 

[Text]  An  aqueous  solution  of  a  synthetic  low-molecular  organic  acid  (NMK)  is  often  used  as  a  cleaning 
agent  during  the  process  of  chemical  cleaning  of  carbonaceous  deposits  on  heat  exchanger  equipment.  Use 
of  this  acid  in  AESs  involves  the  particular  characteristics  of  the  structiiral  materials  in  the  primary 
equipment  ’  .  As  was  shown  in  steps  need  to  be  taken  to  reduce  the  rate  of  corrosion  of  the 
equipment  during  cleaning  in  connection  with  the  aggressiveness  of  the  cleaning  solution  as  compared  to 
the  metal  from  which  the  deposits  are  being  removed. 

One  method  of  reducing  corrosion  in  the  copper  alloys  used  as  the  primary  material  in  the  manufacture  of 
heat  exchangers  is  the  protective  screen  method.  Before  chemical  cleaning,  protectors  —  steel  plates  —  are 
mounted  along  the  loop  being  washed.  According  to  the  data  in  the  literature  one  may  say  that  the  rate 
of  corrosion  of  the  protectors  increases  sharply,  taking  most  of  the  load  of  the  corrosion  process  upon 
themselves,  and  because  of  this  the  rate  of  corrosion  in  the  primary  material  of  the  heat  exchanger,  in 
particular  MNZh-5-1,  should  be  reduced  considerably;  MNZh-5-1  is  used  mainly  in  the  manufacture  of  the 
tubing  systems  for  the  condensers  in  turbines  for  high  power  generating  units. 

However,  we  should  note  that  as  the  steel  protectors  dissolve,  an  additional  significant  amount  of 
trivalent  iron  ions  land  on  the  loop  being  cleaned,  and  the  effect  on  the  nature  of  the  corrosion  process  is 
clear:  there  is  a  dramatic  increase  in  copper  alloy  corrosion,  i.e.  intense  wear  of  the,  heat  exchanger  during 
chemical  cleaning.  Even  if  protective  screening  is  not  used,  steel  structures  alpng  the  loop  —  condenser 
tube  plates,  tubing  intakes  and  drain  pipes  -  may  take  the  place  of  the  protectors. 

As  experience  has  shown  in  the  use  of  periodic  chemical  cleaning,  a  lot  of  cleaning  solution  (and,  hence, 
money)  is  lost  at  the  ends  of  condenser  intakes  and  drain  pipes  and  in  the  places  where  the  tubes  roll 
around  in  the  tube  plates,  i.e.  where  intense  local  corrosional  wear  takes  place. 

With  the  idea  of  reducing  the  impact  of  trivalent  iron  ions  on  the  corrosion  rate  of  copper  alloys  in  acidic 
cleaning  solutions,  we  recommend  the  introduction  of  compounds  that  facilitate  the  transfer  of  trivalent 
iron  to  the  divalent  state  Thiourea  is  widely  employed  for  this  purpose,  in  addition  to  a  number  of 
other  compounds.  However,  the  corrosion  process  in  which  steel  and  MNZh-5-1  alloy  are  present  together 
in  acidic  solutions  with  admixed  reducing  compounds  inhibitory  to  acidic  corrosion  has  not  been  studied 
enough.  There  is  cause  for  interest  in  carrying  out  an  investigation  in  this  area. 
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Figure  1.  A  Graph  of  the  Corrosion  Processes  For  MNZh-5~l  and  Steel  3  Being  Present  Together  in  a  7% 
NMK  Solution,  T  =  333  K:  i  -  MNZh-5-1  corrosion  rate  (No  Steel  3  present);  2-  MNZh-5-1  corrosion 
rate;  3 -effect  of  thiourea  (0.25  g/dm^)  on  the  MNZh-5-1  corrosion  rate;  4  -  6 -  protective  action  of 
acidic  corrosion  inhibitors  (US-1,  R-12,  R-14,  R-15)  on  MNZh-5-1;  7— Steel  3  corrosion  rate  (no  MNZh-5- 
1  present);  3— Steel  3  corrosion  rate;  9  — effect  of  thiourea  (0.25  g/dm^)  on  the  Steel  3  corrosion  rate; 
i3—  i3  —  protective  action  of  acidic  corrosion  inhibitors  (US- 1,  R-12,  R-14,  R-15)  on  Steel  3 

Figure  2.  A  Graph  of  the  Corrosion  Processes  for  MNZh-5-1  and  Steel  3  Being  Present  Together  in  a  7% 
NMK  Solution,  T  =  313  K:  Notation  the  same  as  for  Figure  1. 

The  methodology  for  conducting  a  laboratory  experiment  has  already  been  described  Steel  3  and 
MNZh-5-1  samples  are  simultaneously  placed  in  chemical  beakers  containing  a  7%  NMK  solution,  to 
which  thiourea  is  added  (0.25,  0.5  and  1.0  g/dm^,  tests  3  —  5,  see  Table),  as  well  as  certain  acidic 
corrosion  inhibitors  that  we  synthesized  ourselves  (US-1,  R-12,  R-14,  R-15),  having  a  constant 
concentration  of  0.5  g/dm^  (tests  6  —  9,  see  Table).  In  addition.  Steel  3  and  MNZh-1-5  were  placed  in 
beakers  with  NMK  separately  (test  1)  and  together  (test  2)  without  adding  thiourea  or  inhibitor  in  order 
to  study  the  effect  of  MNZh-5-1  on  the  corrosion  rate  of  Steel  3  and  vice  versa. 

The  Table  shows  the  results  of  the  corrosion  tests  at  298  K,  while  Figs.  1  and  2  show  them  for  313  and 
333  K. 

As  follows  from  the  obtained  data,  the  corrosion  rate  of  Steel  3  in  a  7%  NMK  solution  with  no  additives 
or  MNZh-5-1  is  practically  invariant  over  time.  By  contrast  there  is  even  a  small  decrease  in  the 
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.W  onuTfl 

TEST 
- 1 

Ct.  3 

TIME,h 

% 

A', 

% 

i 

0,09 

0 

0.09 

0 

24 

0,08 

0 

0,08 

0 

48 

0.08 

0 

0,08 

0 

72 

0,09 

0 

0,07 

0 

96 

0.09 

0 

0,07 

0 

120 

2* 

0,13 

0 

0,09 

0 

24 

o.id 

0 

0,08 

0 

72 

0,16 

0 

0,09 

0 

96 

0,16 

0 

0.09 

0 

120 

3 

0,07 

47 

0,17 

-94 

24 

0,10 

26 

0,16 

-91 

48 

0,11 

29 

0,12 

—64 

72 

0,11 

29 

0,14 

-60 

96 

0.15 

8,5 

0,08 

-9 

120 

4 

0,17 

-29 

0,03 

66 

24 

0,15 

-10 

0,02 

87 

48 

0.16 

•  3 

0.02 

83 

72 

— 

— 

96 

0.11 

32 

0,03 

63 

120 

5 

0,12 

9 

0.02 

88 

24 

0,11 

19 

0.02 

87 

48 

0.01 

95 

0.02 

83 

72 

— 

96 

0.06 

63 

0.02 

76 

120 

a 

0,02 

80 

0.01 

89 

24 

0.02 

80 

0,01 

88 

48 

0,03 

79 

0,01 

86 

72 

0.03 

79 

0,01 

92 

96 

0,03 

80 

0.02 

78 

120 

7 

0.02 

83 

0,005 

94 

24 

0,02 

76 

0.006 

93 

48 

0,04 

76 

0,009 

90 

72 

0,03 

78 

0,009 

90 

96 

0.03 

80 

0,007 

92 

120 

8 

0,02 

86 

0,05 

94 

24 

0.02 

86 

48 

0.03 

81 

0.008 

90 

72 

0,03 

78 

0,009 

90 

96 

0,03 

79 

0,009 

90 

120 

9 

0,02 

84 

0,01 

92 

24 

0,02 

86 

48 

0,03 

80 

0,01 

87 

72 

0,03 

78 

0,01 

90 

96 

0,03 

80 

0,01 

86 

120 

10 

0,14 

87 

0,01 

91 

24 

0,06 

60 

48 

0,06 

64 

0,009 

90 

72 

0,06 

64 

0,01 

90 

96 

0,04 

73 

0,01 

89 

120 

*Steel  3  and  MNZh-5-1  are  present  together 


corrosion  rate  of  the  MNZh-5-1  alloy  in  the  same 
situation,  which  agrees  with  the  earlier  work 
The  test  were  conducted  under  static  conditions 
without  stirring  of  the  acid  solution.  The  Steel  3 
and  MNZh-5-1  corrosion  rates  were  found  by 
averaging  five  parallel  tests,  as  is  instructed  in 
the  CEMA  4815-84  standard. 

When  Steel  3  and  MNZh-5-1  are  both  present, 
the  corrosion  process  changes.  First  of  all,  the 
MNZh-5-1  corrosion  rate  “stabilizes”  at  a  higher 
level,  and  secondly,  there  is  observed  a 
substantial  increase  in  the  Steel  3  corrosion  rate 
(by  a  factor  of  1.6  to  2).  In  addition,  there  is  a 
significant  amount  of  reduced  copper  noted  on 
the  steel  plates,  which  is  difficult  to  remove  by 
mechanical  means. 

Introducing  thiourea  into  the  NMK  solution 
turns  out  to  have  a  differing  effect  on  the 
corrosion  process  depending  on  its  concentration. 
For  a  thiourea  concentration  of  1  g/dm^  the 
degree  of  protection  for  Steel  3  (the  decrease  in 
the  corrosion  rate)  fluctuates  between  47  and 
8.5%  by  the  end  of  the  experiment,  while  the 
degree  of  protection  is  chiefly  negative  relative  to 
MNZh-5-1,  i.e.  there  is  an  increase  by  a  factor  of 
1.6  to  1.9  in  the  corrosion  rate.  The  most 
beneficial  thiourea  concentration  in  this  instance 
should  be  recognized  at  0.25  g/dm'^  (test  10, 
Table),  for  in  this  instance  the  Steel  3  corrosion 
rate  was  considerably  reduced.  Nevertheless  we 
need  to  point  out  that  in  all  the  tests  using 
thiourea  (tests  3  —  5,  Table)  there  was  significant 
precipitation  of  elementary  sulfur  and,  in 
addition,  to  some  degree  or  another  all  the  steel 
plates  had  reduced  copper  deposits  on  them. 

The  effectiveness  of  chemical  cleaning  is  reduced 
when  precipitate  falls  into  the  recessed  areas  of 
the  loop,  and  the  presence  of  cupric  deposits  on 
the  steel  surfaces  leads  to  an  increase  in  the 
corrosion  rate  of  the  loop  as  a  whole  and  to  the 
need  for  mechanical  cleaning  or  other  measures. 


The  best  idea,  as  we  have  proposed,  is  to  apply  a  reliable  protective  inhibitor  to  the  entire  loop  to  be 
cleaned^.  It  was  earlier  shown^’®  that  a  number  of  the  inhibitors  we  synthesized  effectively  protect 
against  oxidic  corrosion  in  ferrous  metals  as  well  as  in  the  MNZh-5-1  alloy.  For  reduced  corrosion  rates  in 
the  most  highly  stressed  sections  of  the  loop  in  chemical  cleaning  of  the  heat  exchanger,  it  is  relevant  to 
study  the  inhibitory  process  when  both  Steel  3  and  MNZh-5-1  are  present  in  the  NMK  solution. 
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The  addition  of  0.5  g/dm^  of  some  inhibitor  gets  us  from  a  78  to  an  86%  degree  of  protection  for  Steel  3 
(Table,  tests  6-9)  and  from  83  to  94%  for  MNZh-5-1  alloy.  The  inhibiting  effect  for  higher 
concentrations  of  these  compounds  is  to  be  found  in  We  should  note  that  there  is  a  complete  absence 
of  reduced  copper  on  the  steel  samples,  consequently  the  electrochemical  process  of  reducing  copper  has 
been  completely  suppressed. 

Negligible  concentrations  of  thiourea  lower  somewhat  the  corrosion  rate  of  Steel  3  in  the  7%  NMK 
solution. 

In  order  to  study  possible  synergy,  a  separate  series  of  NMK  solutions  were  prepared  that  had  thiourea 
and  oxidic  corrosion  inhibitors  added  at  the  same  time.  As  was  demonstrated  in  the  investigation  (Table, 
test  10),  in  spite  of  the  absence  of  reduced  copper  on  the  steel  plates,  the  degree  of  protection  of  the 
control  samples  did  not  increase,  i.e.  no  synergy  was  discovered  (no  increase  in  the  inhibitory  effect). 

In  order  to  study  the  effect  of  cleaning  solution  temperature  on  the  corrosion  process,  we  also  performed 
tests  at  the  higher  temperatures  of  313  K  and  333  K,  which  are  sometimes  employed  when  acid  cleaning 
carbonaceous  deposits.  As  follows  from  the  data  (Fig.  1,  2),  the  nature  of  the  corrosion  process  changes 
almost  not  at  all  with  increasing  cleaning  solution  temperature  in  this  case,  i.e.  the  presence  of  thiourea 
facilitates  increased  corrosion  wear  in  MNZh-5-1  alloy. 

Thus  we  have  studied  the  effect  of  both  Steel  3  and  MNZh-5-1  being  present  in  an  NMK  solution  on  the 
corrosion  process.  We  have  remarked  on  the  role  of  protective  screens  in  chemical  cleaning  of 
carbonaceous  deposits  from  heat  exchangers.  The  ineffectiveness  of  using  thiourea  under  actual 
circumstances  has  been  recognized.  The  effectiveness  of  using  a  reliable  protective  inhibitor  has  been 
confirmed. 
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[Text]  A  great  many  studies  have  been  devoted  to  investigating  the  thermal  and  hydraulic  processes 
going  on  in  non-stationary  and  emergency  modes  of  operation  of  VVER-lOOO-equipped  AES.  Most  of  the 
mathematical  models  and  programs  that  have  been  developed  are  directed  toward  reaffirming  the 
operating  safety  of  the  power  generating  units. 

The  problems  that  arise  during  power  generator  operation  place  particular  demands  on  a  mathematical 
model.  We  will  list  a  few  of  these. 

1.  We  need  to  have  a  mathematical  model  of  the  power  generating  unit  as  a  whole,  including  the 
primary  process  systems,  automatic  regulators  and  measuring  devices. 

2.  The  mathematical  model  should  be  able  to  be  readapted  for  each  different  actual  AES  power 
generating  unit  (taking  into  account  the  individual  technical  implementations,  adjustments  to  the 
automatic  regulating  systems,  etc.),  as  well  be  able  to  process  output  data  for  monitoring  by  AES 
technologists. 

3.  The  mathematical  model  should  be  both  fast  and  accurate. 

Satisfying  the  latter  requirement  is  accomplished  by  developing  a  set  of  mathematical  models,  each  of 
which  has  its  limited  area  of  application  while  ensuring  efficiency  and  accuracy.  The  combination  of 
mathematical  rnodels  allows  us  to  have  a  wide  range  of  applicability  of  the  mathematical  model  as  a 
whole  while  maintaining  its  speed. 

The  complex  of  computer  programs  that  we  have  developed  were  designed  on  the  modular  principle.  All 
the  AES  process  equipment  and  physical  processes  were  nominally  divided  into  groups,  each  of  which  was 
described  by  a  separate  functional  module  (subprogram).  Most  of  the  functional  models  had  variant 
implementations  differing  in  their  detail  and  accuracy  of  describing  the  corresponding  object  or  process 
and,  as  a  result,  the  computational  characteristics  of  the  program:  its  speed  and  the  amount  of  computer 
memory  needed.  The  software  package  included  the  following  functional  modules: 

A- a  subprogram  for  neutron  physics  calculations;  Tf- a  subprogram  for  modeling  the  thermal  and 
hydrodynamic  properties  of  the  coolant  in  the  primary  loop  of  the  AES;  TEPL-a  model  describing  the 
temperature  fields  in  the  fuel  elements,  reactor  shell  and  piping;  rP5- joining  the  models  of  the  primary 
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Figure  1.  Structural  Diagram  of  the  KIPR  Software  Complex.  Key:  1  —  Control;  2  Data;  3  Display , 

4  —  Magnetic  Disk;  5  —  Common  Area  in  Computer  Main  Memory 

loop  process  system  and  reactor  core  emergency  cooling  system,  including  the  pumps  for  emergency 
reactor  flooding  and  aftercooling;  PG  — a  model  of  the  steam  generators  and  steam  lines;  — a  model  of 
the  secondary  loop  thermal  equipment,  including  the  turbine,  condenser,  deaerator,  steam  superheater, 
etc.;  and  A  UTO  —  a  functional  subprogram  for  joining  all  automatic  regulating  systems  of  the  primary 
and  secondary  loops  of  the  power  generating  unit.  The  operation  of  the  automatic  systems  is  described 
using  the  real  control  laws  implemented  by  actual  regulators,  their  adjustment  parameters,  dead  zones, 
servomotor  inertia,  etc.  A  A  is  the  model  of  the  AES  power  system. 

To  ensure  operation  of  all  the  program  modules  in  the  context  of  the  overall  mathematical  model,  we 
standardized  the  content  and  form  of  data  exchange  between  program  modules,  and  the  rules  governing 
their  operation.  Ail  data  exchanged  between  the  functional  modules  is  located  in  a  common  area  of  the 
computer’s  main  memory.  A  dedicated  subprogram  monitor  controls  the  operation  of  the  functional 
modules,  which  organizes  calculation  of  the  initial  steady-state  conditions  of  the  power  generating  unit 
and  calls  the  functional  modules  in  cyclic  order  (modeling  a  process  that  is  non-steady  with  respect  to 
time).  The  monitor  also  provides  interactive  operation  with  the  model,  so  that  the  researcher  has  the 
opportunity  to  output  the  model’s  results  in  graphic  and  text  form  to  any  number  of  displays  and  vary 
the  operating  mode  of  the  program  in  an  interactive  manner.  In  addition  to  the  functional  modules  and 
control  subprogram  as  a  part  of  the  software  package,  there  are  also  a  number  of  service  programs, 
including  M/T- a  module  that  accompanies  the  initial  data  needed  by  the  software  package.  Initial  data 
includes  several  thousand  numbers,  listed  in  tables,  laden  with  comments  and  special  indicators.  The 
INIT  service  program  substantially  simplifies  the  preparation  of  initial  data  for  calculations  and  allows  us 
to  construct  a  mathematical  model  for  an  actual  power  generating  unit. 

The  structure  of  control  arid  data  transfer  of  the  software  package  is  shown  in  Fig.  1. 

This  organization  of  the  software  complex  has  the  following  important  advantages  permitting  efficient  use 
of  this  tool  for  operational  tasks: 

•  it  enables  us  to  broaden  the  content  and,  consequently,  to  increase  the  detail  and  accuracy  of  the 
model; 

•  it  gives  us  flexibility  in  the  modeling  systeni  (more  accurate,  slower  models  used  only  in  certain 
circumstances,  or,  for  the  majority  of  calculations,  faster  and  simpler  models); 

•  thete  are  the  service  programs  that  allow  AES  technologists  to  readjust  the  mathematical  model  for 
actual  power  generating  units  and  process  the  output  data  (including  interactively). 

We  will  now  give  a  more  detailed  description  of  the  fundamental  functional  modules  of  the  software 
package.  The  calculational  complex  consists  of  4  interchangeable  versions  of  a  model  of  the  thermal  and 
hydrodynamic  state  of  the  coolant  in  the  primary  loop  of  a  power  generating  unit  (HJ  - 
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1.  The  primary  job  of  functional  module  HI  is  to  model  the  processes  associated  with  large-scale  flows  in 
the  primary  loop  of  the  power  generating  unit  (including  the  Maximum  Design  Accident).  What  is  being 
determined  in  this  case  is  the  correct  description  of  the  effect  of  pressure  wave  propagation. 

One-dimensional  time-dependent  control  of  the  conservation  of  mass,  momentum  and  energy  in  a 
homogeneous  steam-water  flow  is  integrated  in.  It  uses  the  calculational  method  of  which  is  a 
generalization  of  the  Godunov  numerical  scheme^  for  the  case  of  more  complex  (as  compared  to  an  ideal 
gas,  for  which  the  Godunov  method  is  traditionally  employed)  control  of  the  state  of  the  steam-water 
mixture. 

We  note  the  advantages  of  the  numerical  scheme  HI  is  based  on  as  compared  to  other  schemes  for 
hydrodynamic  calculations  (e.g.,  the  Lax-Vendroff  scheme,  or  the  method  of  characteristics):  the  HI 
numerical  scheme  describes  more  accurately  the  wave  flow  of  coolant  for  a  minimal  number  of 
differencing  steps;  conservation  of  mass  and  energy  is  built  in  to  the  numerical  scheme,  which  is  especially 
important  for  investigating  loss-of-seal  in  a  loop;  and  the  numerical  scheme  is  highly  effective  and  stable 
for  supersonic  flows  as  well  as  those  approaching  sound  speed. 

2.  The  primary  job  of  model  H2  is  also  a  description  of  the  processes  associated  with  loss-of-seal  in  the 
primary  loop.  A  one-velocity  two-temperature  model  of  coolant  flow  is  employed.  It  is  assumed  that  the 
temperature  of  the  steam  is  identical  to  the  saturation  temperature  and  it  considers  only  the  most 
significant  effect  -  the  possibility  of  water  making  the  transition  to  a  superheated  state  in  the  initial 
stages  of  boiling. 

The  kinetic  equations  for  steam  generation  include  a  function  for  the  steam  generation  intensity,  which 
follows  from  a  generalization  of  the  large  amounts  of  experimental  data.  It  considers  boiling  at  the  walls 
of  the  tubing  as  well  as  within  the  volume  of  the  superheated  liquid  on  solid  impurities. 

Article  ^  presents  the  calculational  model  on  which  functional  module  H2  is  based  in  more  detail. 

3.  When  modeling  time-dependent  AES  operating  modes  with  characteristic  times  on  the  order  of  tens  of 
seconds,  use  of  mathematical  models  HI  and  H2  is  involves  significant  outlays  of  machine  time.  This  is 
in  connection  with  the  Curant  limit  Ai  <  Aa:/  1 db  c  | ,  which  determines  when  to  ‘‘compress”  the 
propagation  of  acoustic  disturbances. 

We  may  point  out  that  there  is  a  rather  broad  class  of  slow  time-dependent  modes  in  which  wave  effects 
are  not  decisive  and  use  of  mathematical  models  HI  and  H2  is  not  justified. 

Mathematical  model  H3  is  based  on  a  quasi-stationary  calculation  method  for  the  coolant  flow  in  the 
loops  of  a  power  plant  This  method  describes  the  quasi-stationary  flow  class  in  which  variations  in 
velocity  and  pressure  over  time  are  governed  primarily  by  the  time-dependent  boundary  conditions  and 
variations  in  the  external  parameters.  Nondimensional  analysis  shows  that  an  integrated,  simplified 
system  of  hydrodynamic  equations  is  accurate  enough  for  this  class  of  flows:  time-independent  equations 
for  velocity  v  and  pressure  p,  and  time-dependent  equations  for  the  density  p.  The  Curant  stability 
condition  for  this  case  has  the  form  A^  <  Aar/  |  u  |  ,  which  enables  us  to  construct  solutions  with  larger 
time  steps. 

For  description  of  the  effects  of  boiling  and  condensation  caused  by  pressure  changes  over  time,  we  use  a 
constant-volume  idealization,  appropriate  to  which  we  introduce  the  unknown  parameter  which 

characterizes  the  speed  of  pressure  changes  over  the  loop. 
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For  integrating  the  continuity  equations  we  use  a  conservative  numerical  scheme  *  based  on  analytical 
solution  of  the  decay  of  an  arbitrary  density  discontinuity  in  a  two-phase  two-velocity  steam-water  Oow 
The  fundamental  steam-water  flow  regimes  are  considered:  bubble,  slug,  eddy,  disperse,  and  transi  lona 

regimes. 


Thus  functional  module  H3  ensures  conservation  of  mass  and  energy,  takes  into  account  the  expansion  of 
coolant  because  of  changes  in  pressure  over  time,  and  correctly  describes  the  effects  of  velocity 
nonuniformity  with  a  small  number  of  calculation  steps. 


4.  In  the  modeling,  a  fairly  broad  class  of  nonstationary  AES  modes  do  not  require  a  high  deg^e  o 
accuracy  in  descriptions  of  the  coolant  flow  in  the  primary  loop  of  the  power  generating  unit.  The  slow, 
nonstationary  processes  with  characteristic  times  from  hundreds  to  thousands  of  seconds  in  which  the 
coolant  remains  one-phase  are  an  example  of  this. 

We  have  developed  a  simplified  hydrodynamics  mathematical  model  H4  for  description  of  these  regimes, 
and  for  performing  multivariate  quantitative  calculations. 

The  primary  loop  of  a  power  generating  unit  is  represented  by  three  elements  (the  mixing  chambers  and 
the  pressurizer)  with  averaged  parameters.  Changes  in  the  flowrate  of  coolant  over  the  length  of  a  loop 
and  core  channels  are  ignored,  and  we  also  do  not  take  into  account  the  time  it  takes  for  coolant  to  pass 
through  the  loops  and  core  channels.  The  key  to  the  method  is  the  following  assertion: 

dPNKs/*^^  -  ‘^PVKsI^^-  ^ 


From  the  equations  of  conservation  of  mass  and  energy,  we  obtain  a  relation  for  the  rate  of  change  m  the 
pressure  in  the  NKS  and  VKS  [not  further  identified]  with  respect  to  time. 

The  total  system  of  equations  for  determining  the  pressure  and  flowrate  in  the  loop  includes  equations  of 
the  form  Ap  =  (where  4  is  a  known  coefficient)  for  the  loop  and  core  channels,  the  equation  o 
conservation  of  energy  for  the  pressurizer,  and  balancing  equations  of  conservation  of  energy  for  the  loop 
as  a  whole.  We  propose  an  implicit  calculation  scheme  for  solving  these  equations,  which  allows  us  to 
increase  the  time  integration  step  up  to  1  s.  This  makes  the  functional  model  quite  fast  In  additional, 
mathematical  model  H4  provides  enough  accuracy  for  engineering  calculations,  which  follows  troin  a 
comparative  analysis  of  its  operation,  experimental  data,  and  calculations  from  the  models  H3  and  h  . 

The  software  package  includes  two  interchangeable  models  for  neutron  physics  processes  in  the  reactor 
core  (Nos  1  and  2). 

1.  This  model  is  designed  to  describe  the  processes  in  which  spatial  neutron  physics  effects  a.re 
insignificant.  Mathematical  model  No  1  is  based  on  integrating  the  equations  of  point  kinetics  in  a  quasi- 
stationary  approximation  with  delayed  neutrons  being  broken  down  into  six  different  groups.  It  also 
considers  fast  feedback  on  fuel  temperature,  water  temperature  and  boron  concentration.  For  a 
calculation  of  the  reactivity  induced  by  the  control  rods  we  use  an  experimental  dependence  of  the 
integral  effectiveness  of  a  control  rod  bundle  as  a  function  of  the  depth  of  penetration  into  the  core. 

2.  Module  No  2  is  designed  for  a  description  of  the  spatial  distribution  of  the  neutron  flux  in  the  core. 
The  core  is  divided  into  10  calculation  regions,  and  in  each  of  these  the  neutron  flux  is  considered  to  be 
piecewise-continuous.  A  quasi-stationary  one-velocity  diffusion  approximation  is  used,  with  delayed 
neutrons  being  broken  down  into  six  different  groups. 
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medium  is  represented  in  the  form  of  a  piecewise- 
continuous  function  that  is  independent  of  the 
spatial  coordinates  at  the  boundaries  of  each  of  the 
calculational  regions). 

The  very  real  advantage  of  this  model  is  its 
efficiency,  which  makes  it  possible  to  use  neutron 
kinetics  model  No  2  for  multivariate  calculations. 

A  mathematical  model  of  the  horizontal  steam 
generator  has  been  developed  in  two  interchangeable 
variants  (functional  modules  PGi  and  PG2), 
differing  in  the  detail  in  their  description  of  thermal 
and  hydrodynamic  processes,  which  are  realized  in 
steam  generators  in  the  time-dependent  operating 

Figure  3.  Graph  of  Variation  in  Parameters  for 
Switching  Off  One  of  the  Main  Circulating  Pumps  at 
YuUAES:  l  —  Flowrate  of  steam  from  the  “emergency” 
steam  generator;  2  —  Flowrate  of  feedwater  into  it;  3  —  Level 
in  the  steam  generator  with  respect  to  a  baseline  of  1000 
mm;  4  —  difference  in  temperature  between  the  hot  side  and 
cold  side;  5  -  pressure  in  the  steam  generators;  6  -  pressure 

in  the  hot  side: - experiment, - -  — 

calculation  for  the  “emergency”  loop, 

X  X  ““  X  calculation  for  the  “non-emergency”  loop, 

*  ^^vel  readings  without  consideration  of  inertia 


Figure  2.  Graph  of  Pressure  in  the  Second  Loop  (i) 
and  in  the  First  Loop  (2),  and  Variations  in 
Temperature  at  the  Reactor  Inlet  (5)  and  Outlet  (4) 
as  a  Function  of  Time  in  the  Closed  SRK  Turbine 
Mode  (With  Triggering  of  the  Emergency  Shutoff): 

DINAMIKA  program  calculations; 

KIPR  program  calculations 

For  integration  of  the  neutron  kinetics  equations  we 
use  a  method  based  on  expansion  of  the  neutron  flux 
density  function  in  terms  of  eignenfunctions  of  the 
Laplace  operator,  which  allows  us  to  reduce  the 
problem  of  integrating  the  initial  differential 
equations  to  solution  of  a  system  of  linear  algebraic 
equations  in  terms  of  the  average  neutron  flux  over 
the  core  (here  the  breeding  characteristic  curve  of  the 


Figure  3.  Graph  of  Variation  in  Parameters  for  a 
False  Triggering  of  the  Scram  System  at  the  KaAES: 

I  —  Neutron  power;  2  and  3  —  Pressure  in  KD  and  GPK;  4 
and  5  —  temperature  at  the  reactor  inlet  and  outlet, 

— — —  experimental, - calculation 

regimes  of  a  power  generating  unit.  Both  models 
( PGl  and  PGS)  include  a  description  of  the  steam 
generators  and  steam  lines  connecting  them  to  the 
turbines.  We  introduce  a  description  of  the  more 
complicated  model  PG2. 

In  modeling  of  steam  generators,  the  most  complex 
problem  is  determining  the  water  level.  The  PG2 
distributed  model  was  developed  for  a  detailed 
description  of  the  processes  in  the  two-phase  volume 
of  the  steam  generator.  It  takes  into  account  the 
following  primary  effects  that  substantially  influence 
the  thermal  and  hydraulic  state  of  the  steam 
generator  and  the  water  level  readings  in  time- 
dependent  processes:  phase  slipping  in  the  steam- 
water  volume  of  the  steam  generator,  natural 
circulation  of  coolant,  variations  in  pressure  drops 
over  the  submerged  tube  plate,  the  complex 
geometry  of  that  plate,  and  the  presence  of  flow  at 
the  edges. 


The  presence  of  the  submerged  tube  plate  and  natural  circulation  of  coolant  in  the  steam  generator  results 
in  a  number  of  instances  in  the  level  gauge  readings  differing  substantially  from  the  actual  amount  of 
water  in  the  generator.  In  determining  the  level  gauge  readings,  model  PG2  calculates  the  drop  in 
pressure  which  will  be  measured  by  the  gauge.  The  equalizing  vessel  with  pulse  lines  is  also  modeled, 
which  enables  us  to  describe  the  effect  of  needle  inertia  in  the  level  readings.  This  ensures  a  fairly  high 
degree  of  accuracy  in  predicting  level  readings  during  time-dependent  processes. 


The  mathematical  model  for  the  unit’s  automatic  regulation  system  is  implemented  in  only  one  vyiant. 
The  functional  module  A  UTO  models  the  work  of  the  primary  AES  power  generating  unit  regulation 
system  in  normal,  transition  and  emergency  modes,  corresponding  to  complete  or  partial  power 
reductions.  The  model  describes  the  control  system  of  the  primary  loop  of  a  reactor,  the  pressurizer,  the 
secondary  loop  and  turbine.  Analysis  of  the  values  of  controlled  and  observed  variables,  realization 
appropriate  to  the  structure  of  the  regulation  and  control  algorithms,  and  the  formation  and  transfer  of 
vector  control  actions  to  the  other  subprograms  are  all  functions  of  the  A  UTO  module.  The  A  UTO 
module  compares  current  parameters  with  set  values  and  generates  safety  signals  or  interlock  signals,  as 
well  as  diagnostic  information. 


The  modeling  of  the  regulation  system  was  carried  out  in  volume,  allowing  us  to  reproduce  the  operation 
of  the  power  generating  unit  over  the  range  of  powers  and  to  realize  all  normal  and  emergency  modes. 
Operation  of  the  regulators  were  described  by  differential  equations  appropriate  to  the  regulation 
guidelines  that  shape  them.  Static  as  well  as  dynamic  regulator  parameters  were  taken  into  account. 

One  of  the  most  important  stages  in  the  development  of  each  functional  module  included  in  the 
mathematical  model  software  package  is  comparing  the  results  of  its  operation  with  experimental  data 
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and  with  calculations  from  more  accurate  models.  This  kind  of  independent  testing  lets  us  verify  that  the 
approximations  made  in  the  mathematical  models  are  fair,  and  to  refine  the  cross-coefficients  in  the 
correlations  used. 

In  order  to  check  on  the  package  as  an  adequate  mathematical  model  of  the  object,  we  conducted  a 
comparison  of  the  calculations  and  experimental  data  obtained  during  dynamic  tests  in  AES.  A 
comparison  was  also  made  with  the  DINAMIKA  program^. 

Fig  2  shows  a  graph  of  the  changes  in  unit  parameters  when  the  turbine  SRK  are  closed  (with  triggering 
of  the  emergency  shutoffs).  The  graph  in  Fig.  3  illustrates  the  calculational  and  experimental  results  in 
which  one  of  the  four  main  circulation  pumps  that  operate  under  nominal  power  is  switched  off,  which 
happened  at  the  YuUAES  .  The  graph  in  Fig.  4  characterizes  the  results  of  the  false  plunging  of  the 
control  rod  bundle  into  the  core,  as  conducted  at  the  KaAES  No  2  unit. 

The  good  agreement  between  calculation  results  and  dynamic  test  data  (about  30  different  modes  were 
processed  in  all)  opens  up  the  possibility  of  using  the  KIPR  software  package  for  meeting  a  wide  range  of 

needs  arising  in  the  process  of  operating  an  AES,  including  the  development  of  instructional  methods  for 
the  operating  personnel. 
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On  the  Operational  Reliability  of  the  LMZ  VK- 100-90  Turbine 

907F0210C  Moscow  ELEKTRICHESKIYE  STANTSII  in  Russian,  No  2,  Feb  90,  pp  27-30 
[Article  by  Yu.M.  Gofman,  C.Tech.S.,  Sverdiovenergoremont] 

[Text]  VK-100-90  turbines  operate  at  p  =  90  kgf/cm^  T  =  SOO'C  and  at  present  have  logged  about 
300,000  h. 

In  accordance  with  the  technical  reequipment  plan  of  USSR  Minenergo,  all  condenser  turbines  of  this 
brand  are  to  be  written  off,  and  turbines  converted  over  to  central-heating  uses  will  be  modernized,  in 
which  the  replacement  of  the  high-pressure  cylinder  in  the  assembly  is  stipulated. 

This  article  presents  material  about  inspection  of  the  metal  in  the  turbine  as  a  whole,  in  both  the  high- 
pressure  and  low-pressure  cylinders. 

High-Pressure  Cylinder  Housing.  By  inspecting  the  condition  of  the  metal  in  this  turbine  after  200,000  h 
of  operation,  we  have  shown  that  the  high-pressure  cylinder  may  be  operated  reliably  for  up  to  500,000 
hi,  2^  This  conclusion  was  drawn  on  the  basis  of  defectoscope  monitoring,  studies  of  microstructure  zones 
over  the  operating  temperature  regimes  of  the  turbine,  analysis  of  the  fine  structure,  and  study  of  the 
phase  analysis  and  mechanical  properties. 

Recently  we  conducted  additional  studies  with  estimates  of  the  crack  resistance  of  20KhML  steel  and 
investigation  of  the  condition  of  the  turbine  metal  after  275,000  h.  As  is  known,  the  cyclic  ultimate 
strength  and  crack  resistance  criteria  should  be  numbered  along  with  the  static  and  long-term  strength  as 
limiting  criteria  of  the  condition  of  the  turbine  housing. 

The  primary  method  for  estimating  the  crack  resistance  in  the  process  of  operating  the  turbine  is 
magnetic  powder  defectoscopy  (MPD)  at  the  radius  transitions.  After  300,000  h  of  operation  the  housing 
parts  were  subjected  to  monitoring  three  and  four  times.  Thirteen  turbines  were  subjected  to  this 
monitoring,  installed  in  three  different  electric  power  stations.  The  turbines  operated  in  a  base-load 
regime.  The  operating  period  of  all  the  turbines  was  practically  identical.  The  number  of  starts  ranged 
from  200  to  500. 

Table  1  presents  the  data  on  the  number  of  detected  cracks,  and  Table  2  gives  the  depth  characteristics  of 
the  cracks. 

All  the  cracks  occurred  only  on  exterior  surfaces. 
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An  investigation  of  the  cracks  in  the  high-temperature  region  of  the  housing  showed  that  the  cracks  are 
characteristic  of  heat-fatigue,  are  filled  with  oxides,  and  have  pinched,  dead  ends  (Fig.  1).  This  character 
of  the  cracks  testifies  to  their  thermal  origin. 


We  should  note  that  the  damage  in  each  of  the  power  stations  is  not  identical.  Thus  at  one  power  station 
there  are  practically  no  defects  at  all  on  cast  parts,  while  conversely  at  another  there  are  defects  on 
practically  all  the  cast  parts. 

The  most  damaged  area  of  the  housing  is  the  exterior  surface  of  the  chamber  of  the  regulator  stage. 

In  recent  years  VTI  and  TsKTI  have  developed  a  method  for  determining  the  crack  resistance  of  the 
metal  in  cast  turbine  housing  parts  ’  .  The  crack  resistance  characteristics  were  determined  by  this 
method  for  the  metal  of  three  turbine  housings  after  200,000  h  of  operation. 
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The  mechanical  properties  estimated  from  standard  samples  are  given  below  (1  is  for  no  cracks  over  the 
entire  period  of  operation,  2  is  for  an  average  number  of  defects  —  12,  and  3  is  the  maximum  number  of 
crack  formations  observed  over  all  the  13  turbines  -  25). 
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The  hot  hardness  Hg  and  critical  coverage  (Tj^  determined  by  the  VTI  method,  and  the  critical  coverage 
(Tq  determined  by  the  TsKTI  method  are  presented  below. 
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From  an  exa.rriination  of  this  data  it  follows  that  the  crack-resistance  characteristics  for  20KhML  steel  in 
all  cases  exhibit  satisfactory  values. 

A  study  of  the  metal  after  275,000  h  of  operation  at  the  recommended  capacity  shows  that  the  metal  is  in 
satisfactory  condition  .  In  addition  to  the  accepted  methods,  the  high-pressure  cylinder  metal  was 
subjected  to  an  evaluation  of  damage  in  terms  of  the  criteria  proposed  in  The  damage  was 
determined  from  electron  microscope  investigation  of  foils,  electron  fractography  of  low-temperature 
breaks,  and  optical  metallography.  Viewing  a  large  number  of  foil  sections  revealed  that  damage  to  the 
metal  was  low,  and  only  a  small  number  of  individual  micropores  0.1  -0.3  //m  in  size  were  detected  along 
the  boundaries  of  the  ferrite  grains  (Fig.  2). 

We  should  note  that  damage  to  the  metal  as  estimated  from  the  dislocation  structure  of  ferrite  was  not 
taken  into  consideration,  since  the  metal  has  a  large-grain  structure  and  there  is  practically  no  separation 

of  carbides  by  the  main  part  of  a  grain.  Thus  fragmentation  and  formation  of  cell-like  structures  were 
noted  in  the  metal. 

Investigation  reveals  that  the  breaks  are  brittle  and  mixed,  encountering  both  transcrystalline  sections  and 
sections  of  between-grain  breaks.  Only  a  few  fine  pores  were  found  by  examining  the  breaks  on  a 
scanning  microscope,  which  established  that  the  grain  boundaries  are  hardly  damaged  at  all  (Fig.  3). 

Metallographic  analysis  on  an  optical  microscope  detected  no  pores. 

Thus  this  investigation  of  housing  parts  has  shown  that  in  terms  of  all  the  indicators  —  crack  effects,  crack 
resistance,  microstructural  conditions,  mechanical  properties,  phase  composition  and  carbide  formation 
processes,  in  terms  of  damage  the  metal  is  in  satisfactory  condition. 
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High-Pressure  Rotor.  Now  we  will  conduct  a  fairly  detailed  investigation  of  the  structural  condition  of 
the  metal  in  the  rotor  after  extended  operation®,  which  revealed  an  absence  of  substantial  alteration  in  the 
structure  and  properties.  We  followed  the  condition  of  a  large  number  of  high-pressure  rotors  using  non¬ 
destructive  monitoring  methods.  Practically  all  the  checked  rotors  were  free  of  defects.  Defects  were 
detected  in  all  on  only  a  few  rotors.  Investigation  of  the  nature  of  the  defects  showed  that  these  cracks 
were  metallurgical  in  origin  and  located  near  segregated  material;  long-term  operation  does  not  lead  to 
them  opening  any  further. 

Fig.  4  presents  a  typical  crack  4  mm  deep  found  in  the  regulator  stage  area,  which  hardly  grew  at  all  after 
210,000  h  of  operation.  This  leads  us  to  the  conclusion  that  the  level  of  active  stress  in  the  rotor  is 
substantially  less  than  what  has  been  calculated. 

Thus  the  stability  of  the  properties  of  34KhMlA  steel  during  extended  operation  and  the  absence  of 
action  on  the  rotor  of  a  cyclic  stress  nature  that  would  cause  crack  generation  all  point  to  the  possibility 
of  further  reliable  long-term  operation  of  the  high-pressure  rotor. 

Low-Pressure  Housing.  When  following  the  condition  of  housings  after  200,000  h  of  operation,  we  noted 
considerable  corrosional  and  erosional  wear.  Washings  were  found  on  the  cylinder  baffles,  washouts  of  the 
guidevane  roots  along  the  inner  generatrix  and  throat  of  the  diaphragm,  washings  on  the  diaphragm 
joints  and  at  the  guidevane  roots,  and  wear  of  the  flow  separator  (Fig.  5). 

Low-Pressure  Rotor.  Washouts  were  found  on  the  disks  along  the  outside  near  the  rivets  on  the  steam 
intake  side  (Fig.  6). 

Annular  cracks  were  located  on  the  radius  transitions  from  the  circumference  to  the  balance  band.  Radial 
cracks  were  found  near  the  rivet  holes  (Fig.  7).  Cracks  were  found  in  the  key  ways  (Fig.  8).  There  were 
washouts  in  the  hub  of  the  disk  below  the  diaphragm  seals  extending  to  the  disk  apron. 

Thus  this  data  indicates  that  the  operational  reliability  of  the  low-pressure  cylinder  falls  off  sharply  after 
200,000  h  of  operation  due  to  corrosional/erosional  wear  of  the  disk,  diaphragm,  baffles  and  flow 
separator,  as  well  as  corrosional/fatigue  cracking  of  the  disks  in  13  to  19  of  the  stages. 

Conclusions. 

1.  The  high-pressure  cylinder  of  the  VK- 100-90  turbine  may  be  reliably  operated  for  up  to  500,000  h  with 
vigilant  monitoring  after  cracks  begin  forming  in  the  disk  slots  beneath  the  blades. 

2.  The  low-pressure  cylinders  are  not  operationally  reliable  after  250,000  h  of  operation,  requiring  large 
outlays  during  capital  repairs. 

3.  Turbines  that  are  to  be  converted  to  central-heating  uses  do  not  require  replacement  of  the  high- 
pressure  cylinder,  as  has  been  specified  in  the  program,  but  rather  of  the  low-pressure  cylinder. 


Bibliography 

1.  Gofman  Yu.M.,  Kasantseva  N.S.  “Investigation  of  the  Condition  of  Housing  Parts  of  the  K-100-90 
Turbine  After  200,000  h  of  Operation”  in  “Tezisy  dokladov  na  soveshchanii  “Sistema  kontrolya  i 
otsenka  nadezhnosti  i  dolgovechnosti  metall  v  energoustanovkakh”  [Reports  to  the  Conference  on 
“Systems  for  Monitoring  and  Evaluating  the  Reliability  and  Durability  of  Metal  in  Power  Plants], 
Moscow,  VDNKh,  1981 


2T 


2.  Gofman  Yu.N.,  Kasantseva  N.S.  “Reliability  of  VK-100  Turbine  Housing  Parts  After  Long-Term 
Operation”,  ELEKTRICHESKIYE  STANTSII,  No  4,  1984 

3.  Gladshteyn  V.I.,  Sheshenev  M.F.  “On  Reliability  Criteria  of  Cast  Metal  Parts  in  Turbines  Made  of 
ISKhlMlFL”,  TEPLOENERGETIKA,  No  6,  1979 

4.  Rukovodyashchiye  ukazaniya.  Vypusk  44.  Metod  opredeleniya  treshchinostoykosti  materialov 
energooborudovaniya  pri  vysokikh  temperaturakh”  [Guidelines.  No  44.  Method  for  Determining  the 
Crack  Resistance  of  the  Metal  in  Power  Generating  Equipment  at  High  Temperatures],  Leningrad, 
TsKTI,  1981 

5.  “Instruktsiya  po  kontrolyu  za  metallov  kotlov,  turbin  i  truboprovodov.  I  34-70-013-84”  [Instructions 
for  Monitoring  the  Metal  in  Boilers,  Turbines  and  Tubing],  Moscow,  Soyuztekhenergo,  1984 

6.  Gofman  Yu.M.,  Losev  L.Ya.,  Kazantseva  N.S.,  Nevolina  G.S.  “Methods  of  Diagnosing  the  Condition 
of  the  Metal  in  Steam  Piping”,  ENERGETIK,  No  10,  1986 

7.  Gofman  Yu.M.,  Losev  L.Ya.,  “Estimating  the  Degree  of  Damage  in  Metals  Operating  at  Elevated 
Temperatures  Under  Stress”,  MITOM,  No  10,  1985 

8.  Boyko  L.M.,  Novgorodtseva  L.B.,  Filimonova  L.Ye.  “Structure  and  Properties  of  Metal  in  the  Rotor 
of  a  K-100-90  Steam  Turbine  After  Long-Term  Operation”,  ELEKTRICHESKIYE  STANTSII,  No 
11,  1984 


28 


INDUSTRIAL  TECHNOLCXJY,  PLANNING,  PRODUCTIVITY 


UDC  539.319 
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[Article  by  M.  V.  Kvimsadze,  N.  D.  Macharadze,  M.  V.  Ninua,  and  O.  N. 
Chavchadnidze ] 

[Text]  The  situation  frequently  arises  in  design  of  optical 
communication  systems  idien  the  light  guide  changes  from  the 
high-pressure  to  the  low-pressure  zone .  A  sharp  drop  of  ^  external 
pressure,  acting  on  the  fiber,  leads  to  intermittent  variation  of  its 
elastic  characteristics  at  the  interface  of  zones ,  which  in  turn  leads 
to  the  corresponding  variation  of  the  dielectric  constant  of  the  fiber. 
Therefore,  the  emission,  propagated  along  the  fiber,  is  reflected  at  the 
zone  interface.  It  is  physically  obvious  that  energy  losses  of  the 
information  signal,  caused  by  reflection,  will  be  minimal  if  the 
pressure  drop  is  smooth  and  extended  on  some  section  of  the  fiber.  It 
should  be  noted  that  smooth  variation  of  the  stress  tensor  components 
contributes  not  only  to  a  decrease  of  reidiant  energy  losses,  but  also  to 
a  decrease  of  the  probability  that  the  filter  will  chip  on  the  section  of 
the  pressure  drop,  vhich  is  in  itself  an  important  problem. 

To  reduce  the  pressure  jump,  a  transition  region  of  the  fiber  from  the 
direction  of  the  high-pressure  zone  must  be  covered  by  protective 
material,  the  elastic  characteristics  and  configuration  of  which  should 
be  selected  so  that  the  radiant  energy  losses  be  minimal. 

The  distribution  of  elastic  stresses  along  the  fiber  on  the  transition 
region  was  studied  in  this  paper  and  the  conditions  at  which  the 
pressure  .jump  is  minimal  were  determined.  Taking  the  sjTnmetr>  of  the 
problem  into  account,  let  us  write  the  equations  of  elasticity  in  a 
cylindrical  coordinate  system  [1]: 


29 


(1) 


dr 


r  1  ^  i  ^ ' 

1 - (ru,) 

+  1 

\.^  dr 

dz 

^  +  (l_2v) 

dr  dz 


=  0, 


(l~2v) 


1  d 


r  dr 


du. 


dr 


dz 


_L  +  — V) 

r  dr  dz 


=0. 


where  is  Poisson’s  coefficient  of  matter,  u,-  is  the  radial 
displacement  of  the  particle,  and  Uz  is  the  displacement  of  the  particle 
along  the  axis  of  symmetry  of  the  structure,  coinciding  with  axis  z. 

Let  us  discuss  the  boundary  conditions  to  the  equations  of  elasticity: 

a)  It  is  obvious  that  radial  displacements  at  the  fiber-protective 
material  interface  are  identical  for  both  media: 


^Ir  (^u)  (Li)  • 


(2) 


Here  and  further  subscript  "1"  is  ascribed  to  the  fiber,  while  subscript 
2  is  ascribed  to  the  protective  material,  and  ro  is  the  radius  of  the 
f iber . 

b)  Assi^ing  that  the  axlhesion  between  the  fiber  and  protective  material 
is  sufficiently  good,  one  can  assume  that  the  axial  displacements  at 
point  r  =  ro  for  both  materials  are  identical: 


I  ^  ^11*  "■  *  ^  r . 


(3) 


c)  The  radial  stresses  at  the  interface  r  =  rc  in  the  fiber  and 
protective  material  are  identical: 


^}rr  {)  )■  '  ir,,J  . 


(4) 


d)  The  protective  material-external  "internal  pressure"  medium  should 
be  set  equal  at  the  interface  to  the  external  pressure  [2]; 


a2mq«,-=/7„..  (5) 


where  are  the  components  of  the  stress  tensor,  nq  are  the  components 
of  the  external  normal  to  the  surface  of  the  protective  layer,  and  pm 
are  the  components  of  external  pressure. 
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If  the  equation  of  the  generatrix  of  the  cone-like  protective  layer  is 
written  in  the  form  of  r  =  r(z),  boundary  condition  (5)  will  assume  the 
form 


dr 

P  —  ^rz  T— "> 

dz 


i^iz  +  P) 


(6) 


The  following  boundary'  conditions  should  be  added  to  the  conditions  of 
(2)-(4)  and  (6) 


Ur{r—0)  —  0,  li,  (2=0)  =  0. 


(7) 


Let  us  find  the  solution  of  system  (1)  in  the  following  form: 


2®  .  -i 

=  /('■)  +  cp  (r)4-,'?cl)(r).  (8) 

■  <  R 

Substituting  these  expressions  into  equations  ( 1 ) ,  it  is  easy  to 
determine  the  functions  contained  in  formulas  (8): 


/  (ri=c'''>  In  _-f  c<->. 
R 


(p  lr)=c'‘^^ 


F  (/-)  =  - 


R  4(I-v) 
1  R- 


cUi  L  In  -1+  rH) 

R  R  r 


1 


!v  R- 


cW  — £  -£-i _ £ _ Lc<^>-f3  n— v)c<=> 

8  1— V 


3  3— 2v  r-  r  r-i  i  ^  .  «-i 

4  (1-v)  R^  R  R 


O  (r) 
3 


1 


l-2v  8F* 

r3 


4vc<®>— £  ^ _ —  c*''  +  3c*-> 

4  1  —  V 

-(4) 


1  — V 


c">  J:1  In  /rin  L-L]  + 


SF®  R 


4(1— v)  R 


R 


r  R 

-l-r(7)  _  +  c<«>  _ 

R  r 


(9) 


where  r<p  (/=].  2,  /=1,  2  .  8)  the  desired  constants,  determined  from 

the  boundary  conditions,  and  R  is  the  cross-sectional  radius  of  the 
protective  layer  at  z  =  0. 
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Let  us  further  consider  cases  when  Poisson’s  coefficients  of  the  fiber 
and  protective  layer  are  identical  iv\  -  i>2  =  i')  •  Having  used  the  known 
formulas  (see  f2]  for  example)  that  link  (r^q  to  the  components  of  the 
particle  displacement  vector  and  substituting  formulas  (8)  and  (9)  into 
boundary  conditions  (2)-{4),  (7),  and  (8),  we  find  the  following 
relations  for  the  desired  constants: 


c<='>=c<=’>+  ^  p  (4’  +3 VC?)),  c.?^=p  (c'g)  +3  vc<l»), 


c{“)=4')+ 


1— 2v 


(c(i’i  +  3vc?>), 


c?» 


pc?) 


-1-  VC?)  +  ^  (4  vc?>  +  (3-4  V)  3c<?) 


2 - I 


8^2 

1— 2v  2(1— V)  [  R  2) 

1— 2v 


— J-. .  (in  v] 

2(1— V)  [  R  2  l_ 


c?>=c?)+  ^*c?)- 


^  .  c<|)-J_El  c?)  (  In  hS] 
2(1— 2v)  2(l-v)  I,  R  .2) 


(10) 


x>7here  the  following  notation  is  introduced 


Ej—E^ 

Ei+E2(l—2v)  ’ 


(11) 


El  and  Eo  are  Young’s  modulus  of  the  fiber  and  protective  material, 
respectively* 

It  now  remains  to  use  boundary  conditions  (6).  It  follows  from  general 
analysis  of  conditions  (6)  that  a  significant  reduction  of  pressure  on 
the  fiber  at  z  0  may  occur  only  with  a  slow  decrease  of  function  r(z). 
One  can  also  establish  that  conditions  (6)  in  the  region  of  z  ^  R  are 
fulfilled  for  any  z,  if  r(z)  has  the  form 

r-R  exp  I 

The  dimensionless  parameter  a  in  this  expression  determines  the  rate  of 
decrease  of  the  generatrix  of  the  cone-line  protective  layer  and,  as 
follows  from  the  above  arguments,  should  satisfy  the  inequality 
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a  <  1. 


(13) 


Substituting  (12)  into  boundary  conditions  (6)  and  taking  into  account 
relations  (10),  one  can  calculate  the  constants  and  one  can  then 

use  them  to  calculate  the  distribution  along  the  coordinates  of  the 
stress  tensor  components  in  the  fiber.  They  have  the  following  form; 


1 

’’  £.+(1-2.) £7  I-(1-2.)P 


2a(I-2y)p  z*] 
l-(l-2v)§ 


•  h  +  2v^[p+ 


l-(l-2v)p  I 


2a(l  — 2v)p* 
+  l-(l-2v)P 


(14) 


Terms  proportional  to  =Tr  were  not  considered  in  derivation  of  formulas 

(14),  since,  besides  small  parameter  a,  they  contain  small  numerical 
coefficients. 

Analysis  of  formulas  (14)  shows  that  a  significant  decrease  of  pressure 
at  z  -  0  occurs  only  if  the  parameters  of  the  protective  layer  R,  i>  and 
£-2  satisfy  the  inequality 


2(1— v) 


(15) 


Thus,  to  reduce  the  energy  losses  of  the  information  signal  passing 
through  the  optical  system  with  the  section  of  the  pressure  drop,  it  is 
necessary  that  the  parameters  of  the  transition  region  satisfy 
conditions  (13)  and  (15). 
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